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Abstract 

The influence of polydispersity on the magnetization is analyzed in a nonequilibrium situation 
where a cylindrical ferrofluid column is enforced to rotate with constant frequency like a rigid 
body in a homogeneous magnetic field that is applied perpendicular to the cylinder axis. Then, 
the magnetization and the internal magnetic field are not longer parallel to each other and their 
directions differ from that of the applied magnetic field. Experimental results on the transverse 
magnetization component perpendicular to the applied field are compared and analyzed as func- 
tions of rotation frequency and field strength with different polydisperse Debye models that take 
into account the polydispersity in different ways and to a varying degree. 
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I. INTRODUCTION 

The prospect of influencing the rotational dynamics of the nanoscaled magnetic particles 
in a ferrofluid by a macroscopic flow and/or by a magnetic field in order to then observe the 
resulting response via the magnetization and/or via changes in the flow has been stimulating 
many research activities ever since McTague measured |5( the so-called magneto- 

viscous effect. Of particular interest are in this context flows that are shear free on the 
macroscopic scale as in a fluid that is rotating like a rigid body with a rotation frequency, 
say, f2 = Qe z . 

While the colloidal magnetic particles then undergo thermally sustained rotational and 
translational Brownian motion on the microscopic scale they co-rotate in the mean with the 
deterministic macroscopic rigid body flow. However, this mean co-rotation can be hindered 
by magnetic torques on their moments when a magnetic field, say, H = H e x is applied 
perpendicular to the rotation axis e z of the flow. The combination of the externally imposed 
forcing of the particle motion by (i) the rigid body flow in which they are floating and by 
(ii) the magnetic torques on their magnetic moments drives the colloidal suspension out of 
equilibrium. Concerning the magnetic moments, this forcing causes the mean orientation of 
the moments, i.e., of the magnetization M of the ferrofluid to be no longer parallel to the 
internal magnetic field H. Instead, M is pushed out of the direction of H as well as of that of 
H thereby acquiring a nonzero transverse component M y . Here it should be noted that in a 
long cylinder Maxwell's equations imply the vector relation H = Ho— M/2 between the three 
fields when they are stationary and homogeneous but they need not be collinear. However, 
in equilibrium, f2 = 0, the three fields are indeed collinear: the equilibrium magnetization 
M eg (H) = M eg (H)H/H is parallel to the internal field H. 

Recently, the transverse magnetization M y of a slender cylindrical column of ferrofluid 
that was enforced to rotate like a rigid body with constant frequency Qe z in an applied 
homogeneous magnetic field Hoe x was measured as a function of Q and Hq p. These 
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magnetic moments of the larger particles in which the magnetic moments are blocked and 
"frozen" in the particles, i.e., those with effective Brownian relaxation dynamics may be 
rotated by the flow out of the direction of the magnetic field. 

Here we consider poly-disperse models with single-particle Brownian as well as Neel 
relaxation dynamics for the different particle sizes. So we ignore the influence of any dipolar 
magnetic interaction and of any flow induced interaction on the (rotation) dynamics of the 
particles. Thereby collective, collision dominated long-range and long-time hydrodynamic 
relaxation dynamics of the ensemble of magnetic moments are discarded since only the 
individual relaxation of each magnetic moment is considered — albeit in the collectively 
generated internal magnetic field H. 



II. EQUILIBRIUM MAGNETIZATION 

In our experiments we used several ferrofluids out of the APG-series of FerroTec. Their 
saturation magnetization was specified by the manufacturer to be Mf a f = 17507 A/m 
(±10%). This corresponds to a volume concentration <fi » 3.6% of the suspended mag- 
netic material. We have measured the equilibrium magnetization of the ferrofluids with a 
vibrating sample magnetometer (LakeShore 7300 VSM) with a commercial PC user pack- 
age. In order to get information on the particle size distribution of the ferrofluid under 



investigation, we used fits Il3j with a lognormal form of t 



re distribution as well as with 



a regularization procedure |6( based on Tichonovs method w|. Generally the equilibrium 
magnetization M eq (H) as a function of the internal magnetic field H can be approximated 
by a superposition of Langevin-functions 

M°«(H) = J2 W A*AH)]- (2-1) 

i 

Here C(x) = coth(x) — 1/x denotes the Langevin-function that depends on the dimensionless 
Langevin-parameter otj(H) = fiorrijH/kBT and Wj are the so-called magnetic weights, rrij 
refers to the magnetic moment of particles with magnetic diameter dj, i.e., rrij = ~d?M%%l 
with M^ k the bulk-saturation magnetization. From Eq. (j2.1|) we can deduce the initial 
susceptibility Xo — Wjd?- and the saturation magnetization M^ a f = ^ Wj of the 

ferrofluid under investigation. 

Fig. H shows the experimentally determined equilibrium magnetization M eq (H) of APG 
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933 versus internal field H together with fits that were obtained with a lognormal distri- 
bution H and wn, the regulation m e t hod Q. The saturation magnetization of the 
ferrofluid sample was MfJ[ = 19108.6 A/m. From the saturation magnetization the volume 
concentration of the magnetite particles was found to be = Mf a f /M^j = 4.1 %, in rea- 
sonable agreement with the manufacturer's specifications. For the initial susceptibility we 
found the value Xo — 1-09. 

The magnetic weight distributions w(d) resulting from the two fit methods are shown in 
Fig. El 



III. EXPERIMENTAL SETUP 



The experimental setup for measuring the magnetization of a rotating cylindrical column 
of ferrofluid is sketched in Fig. El It is described in more detail in [6] . The ferrofluid is filled 
into a cylindrical plexiglass sample holder with inside radius R = 3.2 mm. This radius is 
so small that for our rotation frequencies the ferrofluid rotates as a rigid body with a flow 
field u(r) = Q x r = Qre^. Here Q is the externally enforced constant rotation rate of the 
sample and e v is the unit vector in azimuthal direction. A homogeneous and temporally 
constant magnetic field Ho = Hq& x is applied perpendicular to the cylinder axis e z . For 
such a combination of enforced rotation and applied field theoretical models allow for a 
spatially and temporally constant nonequilibrium magnetization M that is pushed out of 
the directions of H and H by the flow. 

According to the Maxwell equations the fields H and M within the ferrofluid are related 
to each other via 

H = H - 1m (3.1) 

for our long cylindrical sample and in particular H y = —M y /2 as indicated schematically in 
Fig. El In addition they demand that the magnetic field outside the ferrofluid cylinder 

H- = Ho + i^f2^-M N ) (3.2) 
2 r z \ r r J 

is a superposition of the applied field H and the dipolar contribution from M. This result 
yields a relation between the perpendicular component of the magnetization M y resp. of the 
internal field H y = —M y /2 and the field H y ensor measured by the Hall-sensor outside the 
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sample as indicated in Fig. EJ Considering the finite size of the Hall-sensor, Hy ensor is given 
by 

H , r ^^jy rdx = __£_ Ht . (3 .3) 

In our experimental setup b = 4.75 mm, R = 3.2 mm, and a = 2 mm; here a denotes the 
horizontal extension of the Hall sensor. So, Hy ensor = — 0.386 if^ where H y = —M y /2 is the 
^-component of the internal magnetic field in the ferrofluid. 



IV. MAGNETIZATION DYNAMICS OF A POLY-DISPERSE MODEL 

Comparisons with experimental results showed [6] that theoretical predictions [7] based on 
models 1^,0,0, 3, J^l with a single relaxation time overestimate the magnitude of Hy ensor . 
One reason is that particles with different sizes and different rotational dynamics of their 
magnetic moments contribute differently to the non-equilibrium, flow-induced component 
My of the magnetization and that in particular only the magnetic moments of the larger 
particles with effective Brownian relaxation dynamics may be rotated by the flow out of the 
direction of the magnetic field. 

Therefore, we consider here as a next step poly-disperse models with single-particle Brow- 
nian and Neel relaxation dynamics for the different particle sizes. Such models have been 
used to determine within a linear response analysis the effect of polydispersity on the 
dynamics of a torsional ferrofluid pendulum that was periodically forced close to resonance 
to undergo small amplitude oscillations in a rigid body flow |l3.ll7|. 

We ignore the influence of any dipolar magnetic interaction and of any flow induced inter- 
action on the (rotation) dynamics of the particles. Thereby collective, collision dominated 
long-range and long-time hydrodynamic relaxation dynamics of the ensemble of magnetic 
moments are discarded since only the individual relaxation of each magnetic moment is 
considered — albeit in the collectively generated internal magnetic field H. 

For numerical reasons we use a discrete partition of the particle size distribution. Then, 
without interaction, the magnetization of the resulting mixture of mono-disperse ideal para- 
magnetic gases is given by M = M?, where Mj denotes the magnetization of the particles 
with diameter dj. We assume that each magnetic moment and with it each Mj obeys a simple 
Debye relaxation dynamics that drives them in the absence of flow towards their respective 
equilibrium value M| 9 (H). Then, in the stationary situation resulting from the rigid body 
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rotation with constant ft the Debye relaxation equation for each sub magnetization is given 
by 

flxMj = -[Mj - Mf (H)] . (4.1) 

In the absence of interactions the equilibrium magnetization of each species is determined 
by a Langevin function 

Mf (H) = Xj (H)H = w,C (*?^&h\ § . (4.2) 



Here M^j is the bulk-saturation magnetization of the magnetic material. For the magnetic 
weight Wj(dj) of species j we take the experimentally determined values in the representation 

We should like to draw the attention of the reader to the fact that the magnetization 
equations (|4.1j) for the different particle sizes are coupled via the internal field H: according 
to Maxwell's equations H = H — |M is given in terms of the total M = My. 

In the relaxation rate 1 / Tj we take into account Brownian and Neel relaxation processes 
by adding their rates with equal weight 

111 . , 

T 3 T B T N 

The relaxation times depend on the particle size according to t j b = ^^f (dj + 2s) 3 and 
t n = fo 1 ex P ( 6k bT ^ ■ ^ er ^ ^ S viscosity; s t ne thickness of the nonmagnetic particle layer, 
and K the anisotropy constant. The combined relaxation rate ()4.3|) is dominated by the 
faster of the two processes. Thus, large particles relax in a Brownian manner with relaxation 
times of about some 10~ 3 s, while small particles have the much smaller Neel relaxation times. 
The boundary between Neel and Brownian dominated relaxation as a function of particle 
size d depends sensitively on the anisotropy constant K. This is documented in Fig. 0] for 
the two values K = I5kj/m 3 and K = 50k J/m 3 . For these specific examples the boundaries 
between Neel and Brownian dominated relaxation lie at about d ~ 20nm and d ~ 13nm, 
respectively. 

V. COMPARISON WITH EXPERIMENTS 

For the numerical calculations we take typical values for the ferrofluid APG 933 of Fer- 
roTec that is used among others in the experiments described in j^]: M^j. k = 456k A/ m, 
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7] = 0.5Pa_s, s = 2nm, and fo = 10 9 Hz. Typical values of K lie between I0kj/m 3 and 
50k J/m 3 jig. Il9j. We furthermore use as input the experimental equilibrium magnetization 
M eq (H) of APG 933 shown in Fig. ^and the magnetic weights of Fig. obtained with fits 
to a log-normal distribution or with a regularization method 61. 

n □ 

From the previous work p,|7j we know that single- relaxation time (mono-disperse) models 
predict the maximum of M y resp. of Hy ensor to be located roughly at Qr ~ 1. Furthermore, 
in the experiments j(| done with poly-disperse ferrofluids for frequencies up to Q ~ 3000r ad / s 
mainly the large particles contribute to M y resp. to Hy £nsor since their magnetic moments 
are effectively frozen in the particle's crystal lattice. Only these magnetic moments can 
be pushed out of the direction of the magnetic field by the combined action of thermally 
induced rotary Brownian motion and deterministic macroscopic flow in the rotating cylinder. 
Smaller particles can keep their moment aligned with the magnetic field by the Neel process 
when these particles undergo rotational motion. Finally, the particle diameter that separates 
Neel behavior from Brownian behavior in the size distribution and that thereby determines 
how many particles contribute to M y resp. to the experimental signal H y ensor depends 
sensitively on the anisotropy constant K: The smaller K, the smaller is the number of 
Brownian particles according to Fig. and the smaller is M y resp. Hy ensor . 

The above sketched physical picture is corroborated by Fig. There we compare the 
experimentally obtained H y sensor (Q) (stars) as a function of Q for a representative externally 
applied field Hq = 30k A/ m with various model variants that take into account the poly- 
dispersity to different degrees. This is done for two different anisotropy constants, namely, 
K = 15k J/m 3 and K = 50k J/m 3 as representative examples. However, the three upper- 
most curves refer to single time relaxation approximations, each with r = 2ms 6]: the 
dotted line with crosses is the result of a strictly monodisperse Debye model while the lines 
with diamonds refer to polydisperse models, however, with common Tj = r = 2ms taken in 
Eq. (j4.1)l but different magnetic weights Wj obtained either from a lognormal distribution 
(full line with full diamonds) or from the regularization method (dashed line with open dia- 
monds). The equilibrium magnetization M eq (H) was taken to be the experimental one, the 
distributions were obtained from this experimental M eq (H) by the lognormal ansatz resp. 
the regularization method. This, first of all, shows that models with only one relaxation 
time show roughly the same behavior of M y (Q) irrespective of whether the particle size and 
magnetic moment distributions are polydisperse or not. 
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The set of curves with circles and squares in Fig. |3] refer to truly polydisperse models, 
eqs. ([4.11 - B~3*|) . but different anisotropy constants of the magnetic material [K = 50k J/m 3 
(circles), K = 15k J/m 3 (squares)]. Again, full lines with full symbols were obtained with a 
lognormal distribution while dashed lines with open symbols refer to a distribution resulting 
from the regularization method. Here, one sees that these two distributions with their 
magnetic weights displayed in Fig. El yield very similar results which might not be surprising 
in view of the fact that both seem to reproduce M eq (H) adequately. 

The largest and most important difference between the curves with diamonds (i.e., the 
single-time models) and the curves with circles and squares (i.e., the genuine polydisperse 
models) come from the difference in the anisotropy constants of the magnetic material 
that governs how many particles contribute efficiently as Brownian ones to the transverse 
magnetization M y resp. to Hy ensor : for smaller K the magnetic moments of fewer particles 
being Brownian ones may be rotated out of the direction of the magnetic field by the flow 
in the cylinder. 

The curves for K = 15k J/m 3 yield roughly the same maximal size Hy ensor as the exper- 
iments — they could be fine-tuned even further. But then the location, Q max (Ho), of the 
maxima for different H is still off from the experimental ones as shown in Fig.|Ufa) and[7)^b). 
However, the agreement between the predictions of the polydisperse models of eqs. ()4.1I - 
14. 3|) and the experiments concerning the location Q max (H ) can be improved by allowing 
the relaxation rates Tj of the differently sized particles to depend also on the internal field 
H. To demonstate this, we use for simplicity the form (f| 



with one additional fit parameter 7. Values of about 7 ~ 10~ 4 m/A yield maximum locations 
Q max (Ho) that agree well with the experiments as can be seen in Fig. EJb) and Fig. Efb). 
This generalization of the model ()4.1I - I4.3|) leaves the peak value }j^ ensor [yi max ^ almost 
unchanged, cf, Fig. Ufa). 

However, also this augmented polydispersive model reproduces with fixed values of K 
and 7 the experimental data only in a small range of Q and Ho, cf, Fig. IHland Fig. 
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VI. CONCLUSION 

We have compared the predictions of polydisperse models of the magnetization dynamics 
of ferrofluids with recent experiments measuring the transverse magnetization component 
My of a rotating ferrofTuid cylinder. The models use mixtures of mono-disperse ideal param- 
agnetic gases of differently sized particles. The magnetization dynamics of the models take 
into account the rigid body rotation of the fluid combined with a simple Debye relaxation 
of the magnetic moments of each particle with size dependent Brownian and Neel magnetic 
relaxation times. Thus, in the absence of flow, each magnetic moment and with it each 
sub-magnetization would be driven independently of the others towards its respective mean 
equilibrium value that, however, depends on the internal magnetic field H being collectively 
generated by all magnetic moments. 

The comparison suggests that mainly the large particles contribute to M y since their 
magnetic moments are effectively frozen in the particle's crystal lattice. Only these magnetic 
moments can be pushed effectively out of the direction of the magnetic field by the combined 
action of thermally induced rotary Brownian motion and deterministic macroscopic flow in 
the rotating cylinder. Smaller particles can keep their moment aligned with the magnetic 
field by the Neel process when these particles undergo rotational motion. 

Finally, the particle diameter that separates Neel behavior from Brownian behavior in the 
size distribution and that thereby determines how many particles contribute to M y resp. to 
the experimental signal H y ensor depends quite sensitively on the anisotropy constant K of the 
magnetic material. K determines how many magnetic moments are "frozen" or "blocked" 
in particles and thus can be rotated by the rigid body flow: The smaller K, the smaller is 
the number of Brownian particles with frozen moments, and the smaller is the resulting M y . 
Or, vice versa, a large transverse magnetization can be expected in ferrofluids with large 
anisotropy constants. 

An analysis of the rotation rates £l max (H ) for which M y is largest indicates that the 
agreement between experiments and model predictions can be improved by allowing the 
relaxation rates of the differently sized particles to depend also on the magnetic field H. 
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FIG. 1: Equilibrium magnetization versus internal magnetic field for the ferrofluid AGP 933 of 
FerroTec. Symbols denote experimental data, solid line fit with a lognormal distribution, and 
dashed line fit with a regularization method. See text for further information 
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FIG. 2: Magnetic weights of the 30 here considered particle sizes (d\ = lnm to d^o = 30nm) 
obtained from measurements of M eq (H) in Fig^ by using a lognormal distribution (solid line, 
circles) and by using a regularization method (dashed line, squares). 
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FIG. 3: Schematic plot of the system. The cylindrical sample holder with inner radius R rotates 
with angular velocity ft in the applied static magnetic field Ho perpendicular to ft. The mag- 
netization My is measured with a Hall sensor. M and H denote the magnetization and internal 
magnetic field of the ferrofluid. Both are constant in space and time. 
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FIG. 4: Relaxation times as a function of particle diameter d: Brownian (dot-dashed line), Neel 
(dotted line for K = 15k J /m 3 , dashed line for K = 50/cJ/m 3 ), and the combination (|4.3|) (solid 
line with crosses for K = I5kj/m 3 , solid line with plusses for K = 50k J/m 3 ). 
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FIG. 5: Hy ensor as function of Q for Hq = 30kA/m. The three uppermost curves refer to single- 
time relaxation approximations, each with r = 2ms: monodisperse Debye model (dotted line with 
crosses), polydisperse models with common tj = r and magnetic weights Wj obtained with a 
lognormal distribution (full line with full diamonds) and with a regularization method (dashed line 
with open diamonds). Curves with circles and squares refer to truly polydisperse models [ eqs. (|4.1I 
- 14.3(1 ] with K = 50k J jw? (circles) and K = 15/cJ/m 3 (squares). Full lines with full symbols were 
obtained with a lognormal distribution. Dashed lines with open symbols refer to a distribution 
resulting from the regularization method. In all models the equilibrium magnetization M eq {H) 
was taken to be the experimental one. 
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FIG. 6: Hy ensor as function of O for Hq = 8.6kA/m (squares), 15kA/m (triangles), 30kA/m 
(circles), and 60k A/ m (diamonds). Full symbols denote experimental data. Open symbols refer 
to the polydisperse model [eqs. 1)4.11 - FO)) ] with lognormal distribution and K = lbkJ/m 3 : in 
(a) the relaxation times Tj are independent of H, in (b) they are replaced by t^(H) (|5.1I) with 
ry = W- A m/A. 
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FIG. 7: Maximum magnitude maxHy ensor = H ^ ensor (yi max ) (a) and location of the maximum 
Qjnax ^) ag f unc tions of the external field Hq. Stars show experimental data. Lines refer to the 
results of polydisperse models with a lognormal distribution and K = l5kJ/m 3 : ii-independent 
relaxation times (solid); //-dependencies with 7 = 0.8- 10 _4 to/^4 (dotted), 7 = 10~ 4 m/A (dashed), 
and 7 = 1.1 • 10~ 4 m/A (dot-dashed). In (a) the differences between the lines are too small to be 
seen. 



